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Phase diagrams of binary mixtures of patchy colloids with distinct 
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We investigate the effect of distinct bonding energies on the onset of criticahty of low functionality fluid 
mixtures. We focus on mixtures of particles with two and three patches as this includes the mixture where 
'empty' fluids were originally reported. In addition to the number of patches, the species differ in the type 
of patches or bonding sites. For simplicity, we consider that the patches on each species are identical: one 
species has 3 patches of type A and the other 2 patches of type B. We have found a rich phase behaviour with 
closed miscibility gaps, liquid-liquid demixing and negative azeotropes. Liquid-liquid demixing was found to 
preempt the 'empty' fluid regime, of these mixtures, when the AB bonds are weaker than the AA or BB 
bonds. By contrast, mixtures in this class exhibit 'empty' fluid behaviour when the AB bonds are stronger 
than at least one of the other two. Mixtures with bonding energies tsB = £ab and caa < ^bb-, were found 
to exhibit an unusual negative azeotrope. 



I. INTRODUCTION 

In the last twenty years the study of colloidal phase di- 
agrams, based on spherically symmetric particle interac- 
tions, revealed a rich phenomenology including new crys- 
tal phases, gelation and glass transitionsii^. The phase 
behaviour of binary mixtures is richer and colloidal bi- 
nary mixtures are expected to have a wider range of ap- 
plications in technology as well as in biology. Examples 
include novel candidate photonic crystals, synthesized by 
van Blaaderen and co-workers, through fine control of the 
colloidal charges^ii, and binary mixtures of eye-lens pro- 
teins, studied by Schurtenberger et al, where attractive 
unlike (spherical) interactions were found to be crucial to 
stabilize the mixture, a mechanism relevant in the pre- 
vention of cataract formatioiiri. 

Nowadays colloidal particles can be synthesised in a 
range of shapes and their surfaces may be functionalized 
in a variety of way&§i, with the result that the particle 
interactions become directional or 'patchy'. The prim- 
itive model of patchy colloids consists of hard-spheres 
with / patches on their surfaces. Patchy particles at- 
tract each other if and only if two of their patches over- 
lap. The attraction between particles is short ranged 
and anisotropic: The patches act as bonding sites and 
promote the appearance of well defined clusters, whose 
structure and size distribution depend on the properties 
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of the patches (/ and the bonding energy) and on the 
thermodynamic conditions (density and temperature). 

Sciortino and co-workers established that /, the num- 
ber of patches or bonding sites per particle, is the key 
parameter controlling the location of the liquid-vapour 
critical pointli^. They showed that, for low values of / 
(approaching 2), the phase separation region is drasti- 
cally reduced, and low densities and temperatures can 
be reached without encountering the phase boundary. 
These low density ('empty') phases were shown to be 
network liquids, suggesting that, on cooling, patchy par- 
ticles could assemble into a glassy state of arbitrary low 
density (a gel). 

Remarkably, the results of the simulations of patchy 
colloidal particles are well described by classical liq- 
uid state theories: Wertheim's first order perturbation 
theory^ predicts correctly the equilibrium thermody- 
namic properties; Flory-Stockmayer— theories of poly- 
merization describe quantitatively the size distributions 
of the clusters of patchy particles, including the appear- 
ance of network (percolated) fluids. 

In subsequent work we addressed, explicitly, the inter- 
play of the entropy of mixing and the entropy of bonding 
in the phase behaviour of models of binary mixtures of 
patchy particles. We focused on mixtures of particles 
which differ in the number of patches or functionality 
only. We found that, within this class of mixtures, the 
difference between the functionality of the particles is the 
key parameter controlling the phase equilibria of the sys- 
tem. In particular, if one species has more than twice the 
number of bonding sites of the other, a phase transition 
between two network fluids appears and the topology of 
the phase diagram may change. The miscibility at high 



2 



pressures is also controlled by this difference and closed 
miscibility gaps are always present when the difference, 
between the functionality of the particles, is greater than 
oneii. The conditions required to realize 'empty' phases 
were shown to differ from those proposed by Sciortino 
and co-workers, and no 'empty' liquids were found when 
the functionality of one of the species exceeds 4ii. 

The primitive model of patchy particles was general- 
ized to include different types of patches (i.e. more than 
one bonding energy). This was achievedi^iiS. by consid- 
ering a single component fluid where the particles have 
two A sites (with bonding energy caa), and one B site 
(with bonding energy cbb)- Unlike sites also interact 
with bonding energy cab- This model allows a deeper 
understanding of the onset of criticality in low function- 
ality systems: The detailed fashion in which the critical 
temperature vanishes, as the bonding energies decrease 
towards zero, depends on the order in which the limits 
are taken, which in turn determines the type of network 
that is formed. 'Empty' phases are low-density struc- 
tured (percolated) fluids and the emergence of criticality 
is related to this structure. Thus, the generalized model 
provides a reference system for the microscopic descrip- 
tion of the competition between condensation and self- 
assembly of equilibrium structured fluids as the theoret- 
ical predictions based on Wertheim's theory have been 
confirmed, recently, by Monte Carlo simulation a^^i^^ . 

In the present work, we proceed to investigate the effect 
of distinct bonding energies on the onset of criticality of 
low functionality fluid mixtures. We focus on mixtures of 
particles with two and three patches as this is the type of 
mixture where 'empty' fluids were originally reported. In 
addition to the number of patches the species may differ 
in the type of patches or bonding sites. For simplicity, 
we assume that the patches on each species are identical, 
i.e. we consider mixtures of particles with two B patches 
and particles with three A patches: 2b — i a mixtures. 

As for pure fluids, we will use Wertheim's ther- 
modynamic perturbation theory (in its extension to 
mixtures)^ and Flory-Stockmayer's theory of percola- 
tion (generalized to mixturesii). The colloids are mod- 
elled by equisized hard spheres and we vary the strength 
of the bonding energies caa, ^bb and tAB (Fig- tt])- We 
recover the 2a — 3a mixture investigated previousljiii 
when all bonding energies are equal. This is a type I mix- 
ture that is completely miscible above the critical pres- 
sure of the less volatile component. The 'empty' fluid 
regime is approached monotonically: The critical pack- 
ing fraction and temperature decrease towards zero as the 
fraction of 2a particles approaches one {i.e., the average 
functionality decreases towards two). 

We found a rich phase behaviour with closed miscibil- 
ity gaps, liquid-liquid demixing, negative azeotropes and, 
from a topological point of view, three different types of 
mixtures: type I, I-A and Illii. Liquid-liquid demixing 
is the mechanism which prevents this class of mixtures 
from exhibiting 'empty' fluid behaviour and this occurs 
when the AB bonds are weaker than the AA or the BB 




FIG. 1. Schematic representation of the mixture under con- 
sideration illustrating the different bonding energies. 



bonds. By contrast, the mixtures exhibit 'empty' fluid 
behaviour when the AB bonds are stronger than at least 
one of the other two. Furthermore, mixtures with bond- 
ing energies ebb — ^ab- and eaa < ^bBi exhibit an 
unusual (negative) azeotropic behaviour. 

We found an exception where the AB bonds are 
stronger but the mixture fails to exhibit 'empty' fluid 
behaviour: the mixture with cab = ^aa and sbb ~ 0. 
Near the critical point of this mixture there is demixing 
and two network fluids with finite packing fraction coex- 
ist, as the pressure vanishes. This behaviour resembles 
that found previouslyii for 2a — 5a mixtures where the 
entropy of bonding drives a demixing transition, which 
preempts the 'empty' fluid liquid- vapor phase transition. 

The remainder of the paper is organised as follows. In 
section In] we present the model, summarise Wertheim's 
theory for binary mixtures (Sec. Ill Ap and the Flory- 
Stockmayer theory of percolation for binary mixtures of 
patchy particles (Sec. IIIB| . In section ITTll we present 
the results: phase diagrams (including percolation lines) 
and critical properties of several representative mixtures. 
Finally, in Sec. IIVI we summarize our conclusions and 
suggest lines for future research. 



II. MODEL AND THEORY 

We consider a binary mixture of A^i and A'2 equisized 
hard spheres (HSs) with diameter a. Particles of species 
1 are decorated with 2 patches of type B, and particles 
of species 2 with 3 patches of type A: a 2^ — 3a binary 
mixture. The bonding sites are distributed on the parti- 
cle surfaces in such a manner that two particles can form 
only one single bond, involving two bonding sites only, 
one in each particle. In addition, there is a minimum dis- 
tance between the bonding sites to ensure that no sites 
are shaded by nearby bonds. 



A. Helmholtz free energy: Wertheim's thermodynamic 
perturbation theory 

A description of Wertheim's perturbation theory for 
pure fluids and mixtures can be found elsewhere^i^!^. Here 
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we summarise the theory for the mixture under consider- 
ation (the theory for binary mixtures with an arbitrary 
number of patches of different types is discussed in^) . 

The Helmhohz free energy per particle, fn, is the sum 
of contributions from a reference system, a mixture of 
hard spheres {/hs), and a perturbation due to the bond- 
ing interactions {ft)- 

fH = Fh/N = fHS + fb, (1) 

where N = Ni + N2 is the total number of particles. The 
HSs free energy may be written as the sum of ideal-gas 
and excess contributions: fns ~ fid + fex- The ideal-gas 
free energy is given by 

= Inr? - 1 + ^ ln(a;(^)VO, (2) 

1=1,2 




FIG. 2. Schematic representation of a tree-like cluster in a 
2b — 3a binary mixture. Pa^~, is the probability of bonding 
a site a to a site 7. 



with P = kT the inverse thermal energy, Vi the thermal 
volume and x''-' = Ni/N the number fraction of each 
species, r] = Vgp is the total packing fraction, Vg = 7r/6cr^ 
the volume of a single particle and p the total number 
density. The excess part, which includes the effect of 
the excluded volume, is approximated by the Carnahan- 
Starling equation of statsi^ (note that both species have 
the same diameter): 



fife. 



4:1] - Srf 



(3) 



The bonding free energy has two contributions: the 
bonding energy and an entropic term related to the num- 
ber of ways of bonding two particles. In the framework of 
Wertheim's thermodynamic first-order perturbation the- 
ory it is given by 



■11.16 



:[21nXs-XB + l]+(l-a 



(4) 

where x = x^^^ is the composition of the mixture (cc'^' = 
1 — x) and Xa is the probability that a site of type 
a = A, B is not bonded. The latter are related to the 
thermodynamic quantities through the law of mass ac- 
tion: 

Xa = 1- 2,jxXbXaAab - 3?7(1 - x)XIAaa 
Xb^I- 2t]xXIAbb - 3r?(l - x)XaXbAab (5) 

The set of parameters A^^ characterise the bonds be- 
tween sites a and 7. For simplicity, the interaction be- 
tween bonding sites is modelled by square well potentials 
with depths e^-y. Assuming that all bonds have the same 
volume, Vb, and using the ideal-gas approximation for 
the pair correlation function of the reference HS fluid, 
we findii 



Act = — [expipeaj) - 1] . 



(6) 



three bonding energies, eAA, ^ab and cbb, are equal and 
vary the third. The two identical energies set the energy 
scale, e. In what follows we use scaled-bonding energies 

We use the Gibbs free energy per particle (17 = p/ p + 
fn, where p is the pressure) to obtain the equilibrium 
properties of the mixture. At fixed values of the compo- 
sition X, pressure p, and temperature T, g is minimised 
with respect to the total density p, subject to the con- 
straints imposed by the law mass action. We use a stan- 
dard Newton-Raphson method to minimise g, and solve 
the law mass action simultaneously by a Powell hybrid 
method. 

A standard common tangent construction on g{x) is 
used to determine the coexistence points (which is equiv- 
alent to imposing the equality of the chemical potentials 
in the coexisting phases) while mechanical and thermal 
equilibria are satisfied by fixing the pressure and the tem- 
perature, respectively. 

Finally, the critical points are computed numerically 
by determining the states which satisfy the spinodal con- 
dition and the vanishing of the third-order derivative in 
the direction of largest growtbi^. 



B. Theory of percolation 

Recentljiii we generalised an extension of the Flory- 
Stockmayer random-bond percolation theorjsiS proposed 
by Tavares et. alt^ to binary mixtures with an arbitrary 
number of distinct bonding sites. Consider a tree-like 
cluster {i.e. with no loops) as schematically illustrated 
in Fig. [21 The number of bonded sites a at the level 
2 + 1 is related to the number of all types of bonded sites 
in the previous level through the recursion relations (see 
referencesiiii^ for details) : 



As in previous worksi^, we set the bond volume to = 
0.000332285(7'^. We consider mixtures where two of the 



rii+i^A = 2pA^Ani,A +PB^Ani^B 
rii+i.B = 2pA^Bni,A + Pb^bit-lb, (7) 
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where Pa^j is the probability of bonding a site a to a 
site 7. Then, the probabihty of finding a bonded site a 
is 

Pa = PA^A + PA^B 

Pb = Pb^a + Pb^b , (8) 

which may be related to the thermodynamic quantities 
through the law mass action, since 

Pa^l-X^, a = A,B. (9) 

A term-by-term analysis of Eqs. ([S]) and ([5]) gives 

Pa^a = (1 - x)3r]X^AAA 

PA^B = x2r]XAXBAAB 

Pb^b = x2r]XgABB 
PB^A = (1 - x)3t]XaXbAab- (10) 

In order to find whether the system is percolated or 
not, we write the equations ([7]) in matrix form 

fii = f'no, (11) 

where fij is a vector with components n^.Q and T is a 
2x2 square matrix with entries 

f = (f^^^ P^^A, (12) 

\2pA^B PB^B J ^ ' 

which may be diagonalized. The progressions defined by 
Eq. (jlip converge to if the largest (absolute value) of 
the eigenvalues of T is less than unity. In other words, 
percolation occurs when |A±| = 1 for any of the two 
eigenvalues A± of T. 

III. RESULTS 

Before describing the results for 2a — 3a mixtures we 
set the graphical code/s used in the figures. Phase di- 
agrams are presented, in the temperature-composition 
plane at constant pressure, in Figs. [3l El and [H In all 
the figures percolation lines are plotted as solid red lines: 
Below the percolation line the fluid is percolated in the 
sense that there is a non-zero probability of finding an 
infinite cluster. We call these networks fluids. The liquid 
side/s of the binodal lines is/are always percolated; bin- 
odal lines are depicted as black solid lines; shaded areas 
are two-phase regions; empty circles are critical points; 
black squares are azeotropic points. 

The critical properties of the mixtures are analysed by 
means of pressure-temperature projections (Figs. 21 13 
and llOp and critical temperature vs critical packing frac- 
tion representations (Figs. O [51 and [TT|) . In all cases: 
Solid lines depict the liquid-vapour transition and empty 
circles the critical point of the pure fluid (3^ fluid). 
Dashed and dotted lines are the critical points of the mix- 
ture: When the mixture exhibits 'empty' fluid behaviour, 
at low pressure, the critical line is dashed, otherwise it is 
shown dotted. 



A. 2a — 3a mixture 

We start by reviewing the results for the mixture with 
identical bonding interactions, which was investigated by 
grand-canonical Monte Carlo simulations^ and using the 
current theorj*ii. At low pressures, the pure 3a fluid un- 
dergoes LV liquid-vapour condensation, below the criti- 
cal temperature. The transition involves two fluids with 
different densities and fraction of unbonded sites. By 
contrast, there is no LV transition in the pure 2a fluid, 
as the particles with two bonding sites can form linear 
chains only and the absence of branching prevents the 
fluids from condensing. For the mixture, the simulation 
and theoretical results indicate that the critical packing 
fraction decreases continuously towards zero as the pres- 
sure vanishes, exhibiting 'empty' liquid behaviour, that 
is, there is a fluid phase with arbitrary low packing frac- 
tion. 

The temperature-composition phase diagrams are il- 
lustrated in Fig. 121 panels (61) and (62) at two distinct 

(3) 

pressures below pc , the critical pressure of the pure 
3 A fluid (the mixture is completely miscible above Pc )■ 
When 2a particles are added the LV phase transition 
shifts from x ~ {3a fluid) to finite values of the compo- 
sition. The two-phase region, which is always bounded by 
a lower critical point, increases as the pressure decreases. 
Near the a; = axis the slope of the binodal is negative, 
indicating a decrease in the stability of the liquid phase 
as the fraction of 2a particles increases. The addition 
of 2a particles reduces the probability of branching, the 
mechanism which drives condensation, and consequently 
the region of stability of the liquid phase decreases. 

The pT projection of the critical line (see Fig. \^ starts 
at the critical point of the less volatile fluid (3a fluid) and 
decreases monotonically to p —> and T 0. Topolog- 
ically this is a limiting case of type I mixtures, with no 
liquid-liquid demixing, according to the classification of 
van Konynenburg and Scottii), where one of the fluids 
does not undergo LV condensation (2a fluid). 

In what follows, we consider generic 2b — 3a mixtures 
{i.e. where species 1 is decorated with 2 patches of type 
B and species 2 with 3 patches of type A) which are 
characterised by three bonding interactions: e*AA- ^*bb 
and e*AB- We proceed by setting, in turn, two of them 
equal and vary the third. 



B. 2s - 3a mixtures: eaa = eas 

We start by setting e^^ = e^^ = 1. Tx phase di- 
agrams are depicted in Fig. O at two different pres- 
sures: pI = piVs/e = 1.05 X 10^"' (left column) and 
P2 = 1-05 X 10~^ (right column), both below pi'^\ the 
critical pressure of the pure 3a fluid. pT projections are 
plotted in Fig |31 and critical properties are illustrated in 

Fig. m 
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FIG. 3. Phase diagrams of mixtures of particles with 2 patches of type B (species 1) and 3 patches of type A (species 2) in the 



scaled-temperature vs composition {x 



plane at constant pressure. The bonding interactions are: e^^ = e\g — 1 and: 



row 1 (ai and 02) t'gg = 1.1, (61) and (62) egg — 1, (ci) and (C2) e*gg = 0.5, (di) and (6(2) e*gg — 0. Left column: Pressure 
pi = piVs/e = 1.05 X 10~*, right column: Pressure p2 ~ P2Vs/e = 1.05 x 10~^. 



1- ^BB > li ^*AA — ^AB — 1 

When e|jg > 1 there is liquid-liquid demixing at low 
temperature for any value of the pressure. The demixing 
involves two structured (percolated) fluids: One rich in 
3a particles, and the other rich in 2b particles. The phase 



separation is driven by the presence of AB bonds as they 
reduce the number of BB bonds which are energetically 
favourable. The demixing at low temperatures preempts 
the 'empty' fluid regime. 



We chose e 



BB 



1.1 as a representative example of this 
class of mixtures. At pressures higher than pi the Tx 
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FIG. 4. Pressure-temperature projections of the phase dia- 
grams of binary mixtures with e\j^ = e\g — 1. The sohd 
curve is the hquid- vapour transition of the pure 3a fluid, 
which ends in a critical point (open circle). Dashed and dot- 
ted lines are the critical lines of the mixture. 

phase diagram (not shown) consists in a demixing region 
bounded above by an upper critical point. Sightly below 

(3) 

Pc there are two regions of phase separation bounded by 
critical points (panel (ai) of Fig. LL demixing at low 
temperatures and LV at high temperatures. By reducing 
the pressure below a certain value {j)/p^c'^ ~ 0.27) the two 
regions merge giving rise to a very large demixing region 
without critical points (panel (02) of Fig. [S]). 

The pure 2b fluid does not undergo LV condensation, 
and thus the phase diagram of these mixtures does not 
fit exactly the standard classification of van Konynen- 
burg and ScotlJi. Nevertheless, it may be understood as 
a limiting case of type III mixtures (see the pT projection 
in Fig. 2]). The critical line, which changes its character 
continuously from LV to LL, starts at the critical point 
of the less volatile fluid, exhibits a minimum, and con- 
tinues with negative slope to high pressures. This line 
is expected to intersect the solid phase boundaries that 
were not considered here. The topology of the phase dia- 
gram remains the same for other mixtures with e|j^ > 1 , 
but changes in the critical line, such as the absence of 
the minimum, may occur. 

The critical temperature and packing fraction are plot- 
ted in Fig. [51 As the pressure increases the packing 
fraction increases while the critical temperature and con- 
centration tend asymptotically to kTc/e — > 0.037 and 
Xc 0.757. 

■ Ess — u, e^A — eas — ^ 

Let us continue by considering the extreme case where 
the B patches do not interact: e*^^ = 0. The mixture 
is completely miscible at pressures above pc . Slightly 
below Pc the Tx phase diagrams (an example is plotted 



in panel (di) of Fig. [3]) are similar to those of the refer- 
ence mixture {2a — 3a): The LV phase transition shifts 
from .T = {3a fluid) to finite values of the composition. 
The differences appear at lower temperatures. As a re- 
sult of the absence of bonding between 2b particles, there 
is no condensation when the fraction of 2^ particles in 
the mixture is sufficiently large. When the temperature 
decreases, at any value of the pressure, the percolation 
line tends asymptotically to cc — > 0.75 (with Xa — > 
and Xb 0.5) rather than a; — >■ 1 (with Xa and 
Xb 0) as in the reference mixture. Differences in the 
phase diagrams are noticeable as the pressure is further 
reduced (see panel (^2) in Fig. [3]). Near x = 0, there 
is a LV transition similar to the transition of the ref- 
erence mixture, but as the temperature is lowered, the 
two-phase region changes drastically and ends in a lower 
critical point. Topologically, this is still a limiting case 
of type I mixtures {pT projection is depicted in Fig. 21). 
The critical temperature approaches zero when the pres- 
sure vanishes, as in the reference mixture, but the criti- 
cal packing fraction tends asymptotically to r]c = 0.0091 
(see Fig. and Xc — > 0.699 rather than ijc ^ and 
Xc — > 1. Thus, there is no 'empty' liquid regime when 
the B patches do not interact. The difference between 
the AA and BB bonding energies drives demixing be- 
tween two network fluids, which preempts the LV phase 
transition at low pressures. 



3. < e*BB < 1, ^AA = ^*AB = 1 

These mixtures retain features of the previous two 
classes, those with e*BB = 1 and those with e*gg = 0. 
An example of the temperature-composition phase dia- 
grams at constant pressure is shown in the third row of 
Fig. Ofor eBB = 0-5, examples of the pT projections are 
depicted in Fig. |4l and the critical properties are illus- 
trated in Fig. |5l 

At intermediate temperatures, the BB interactions are 
frozen (in the sense that the probability of BB bond- 
ing is much smaller than the probability of AB or AA 
bonding) and the mixture behaves as the mixture with 
eBB = 0. Eventually, at low enough temperatures, all 
bonds are established and the mixture behaves as the 
reference 2a — Sa mixture. In particular, the 'empty' liq- 
uid regime {Tc — > and ?7c — > as the pressure vanishes) 
is observed. In all cases these mixtures are limiting cases 
of type I phase equilibria where only one component un- 
dergoes LV condensation. 



C. 2s — 3a mixtures: eas = ess 

We continue the analysis by setting e\g = e^g = 1 
and varying the bonding energy of AA bonds. 
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case of type III phase equilibria (see the pT projection 
in Fig. [7] (a)) where only one fluid undergoes LV con- 
densation. The critical line starts at the critical point 
of the pure 3^ fluid, exhibits a pressure maximum fol- 
lowed by a pressure minimum at lower temperatures and 
it continues with negative slope to high pressures, where 
it is expected to intersect the solid phase lines. The crit- 
ical density increases monotonically as the pressure is 
reduced (a representation of Tc vs rjc is shown in Fig. 
[5] (a)). Topological changes are not expected for other 
values of > 1, although minor modifications of the 
critical line may occur. Thus, there is no 'empty' fluid 
regime in this class of mixtures, as it is preempted by LL 
demixing at low pressures. 



FIG. 5. Scaled-critical temperature vs critical packing frac- 
tion of mixtures with = e^g — 1. The open circle denotes 
the critical point of the pure 3a fluid. The inset is a zoom of 
the case egg = 1.1 where rjc ^ 1. 



■ ^AA > ^1 ^AB — ^BB — ^ 

We selected e^^ = 1.1 to illustrate the behaviour of 
this class of mixtures. Temperature-composition phase 
diagrams at different pressures are plotted in Fig. [n](ai) 
pvs/e = 5.24x 10~^ (well above the critical pressure of the 

pure substance Pc'^''), (02) pvs/e — 3.142 x 10~^ (sightly 

above pi^^), and (03) pi = pVs/e = 1.05 x 10^"* (sightly 
(3) 

below Pc )■ The behaviour is similar to the mixture with 
e*gg = 1.1 and e\g = e^^ = 1, analysed in Sec. IIIIB II 
In addition to the LV condensation, there is a LL demix- 
ing region at low temperatures for any value of pressure. 
Well above pc the LL demixing region is bounded by an 
upper critical point. Fig. [6] (ai). As pc is approached 
from above, the LV condensation appears in the form of 
closed loops of immiscibility (02). By further reducing 
the pressure, the two-phase regions grow and eventually 
merge giving rise to a very large demixing region, which 
intersects the x — axis at the temperature of the LV 
transition of the pure component (03). The percolation 
analysis shows that liquid phases are always network flu- 
ids while the vapour phase may or may not be percolated. 
Closed loops of immiscibility, when present, are entirely 
percolated. 

The origin of the LL demixing at low temperatures is 
the same as that observed in the mixtures with — 1.1 
and = e\j^ = 1. As the temperature is reduced the 
fractions of unbonded sites, Xa and Xb, decrease and 
there is demixing between a phase where the AA bonds 
are predominant (small values of composition) and an- 
other where the BB bonds predominate (large values of 
composition). AB bonds are penalised as they reduce 
the probability of A A bonding, which is energetically 
favourable. 

This mixture is topologically classified as a limiting 



2. < e*AA < 1, Kab = (*BB = 1 

In this class of mixtures it is entropically and ener- 
getically favourable to form of AB bonds and therefore 
these mixtures are expected to exhibit 'empty' fluid be- 
haviour. Let us start by reducing the interaction between 
A patches to = 0.85. Temperature-composition 

phase diagrams at constant pressure are depicted in the 
second row of Fig. |6l The mixture is completely miscible 
at pressures above pi ■ By reducing the pressure below 

(3) 

Pc we find that the LV condensation of the pure 3a fluid 
at the a; = axis is shifted to flnite values of the compo- 
sition and ends in an upper critical point. Panel (61) of 
Fig. [S]is an example. By contrast to the mixtures investi- 
gated previously, the slope of the binodal in the Tx plane 
is positive, indicating an increase in the stability of the 
liquid phase when 2b particles are added to the mixture. 
The addition of 2b particles reduces the probability of 
branching (or equivalently the entropy of bonding) and 
a negative slope of the binodal, signals a decrease in the 
stability of the liquid phase. Here, however, the loss in 
the entropy of bonding is overcome by the gain in bond- 
ing energy as the AB bonds are energetically favourable. 
Note that mixtures where < e|j^ < 1, = = 1 

(analysed in Sec. IIIIB 3[) are similar but not exactly the 
same. In that case the bonding energy of the mixture 
is also minimised by the formation of AB bonds, with 
respect to BB bonds, but it is higher than the bond- 
ing energy of a pure 3^ fluid, whereas in the mixture 
analysed in this section (e^^j < 1, = e|j^ = 1) it 

is possible to achieve mixed states with bonding energy 
lower than that of the pure 2>a fluid. 

As pressure is further reduced, a negative azeotrope ap- 
pears (see panel (62) of Fig. [H]). The second LV branch 
leaving the azeotrope at high composition grows rapidly 
as the pressure is reduced (panel (63) in Fig. [5]). This LV 
branch is bounded by a critical point with critical den- 
sity and critical temperature that decrease continuously 
towards zero as the pressure vanishes (see Fig. |5] (b)). 
Therefore, as expected, the mixture exhibits 'empty' fluid 
behaviour. 

Topologically this is again a limiting case of type I 



8 



X 

0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 1 




FIG. 6. Phase diagrams in the scaled-temperature vs composition (x — ') plane at constant pressure of 2b — 3a binary 
mixtures with — = 1. First row; — 1.1, pressures; p* = pvs/e = 5.24 x 10^^ P* ~ 3-14 x 10~* {0-2), 

p* =pl^ 1.05 X 10"" (as). Second row; e^A = 0.85, pressures; p* = 5.24 x 10"* (61), p* = 3.14 x 10"* (62), p* = 1.31 x 10"* 
(63). Third row; e^aO.S, pressures; p* = 4.19 x 10"* (ci), p* = 2.83 x 10"* (C2), p* = 2.09 x 10"* (ca). Fourth row; e^A = 0, 
pressures: p* = 2.88 x 10~* (di), p* = 2.09 x 10"* (^2), P* = P2 = 1-05 x 10"* (ds). 
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FIG. 7. Pressure-temperature projections of the phase dia- 
grams of different binary mixtures with e\g — e*gg = 1. 
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FIG. 8. Scaled-critical temperature versus critical packing 



mixtures with the addition of a negative azeotrope, usu- 
ally referred to as type I-A phase equilibriaiii (see the pT 
projection in panel (b) of Fig. [7]). 

We proceed by reducing the AA bonding energy to 
^AA ~ • Three Tx phase diagrams at representative 
pressures are depicted in the third row of Fig. [S) The 

(3) (3) 

mixture is completely miscible at p > pc . Below pc 
there is a range of pressures in which the phase diagram 
consists only of the LV phase transition, shifted from 
the X = axis (panel (ci) in Fig. [6]). As before, the 
slope of the binodal is positive. There is a remarkable 
maximum in the percolation line at temperatures higher 
than the LV condensation of the pure component. Thus, 
by increasing the number of particles with two patches it 
is possible to find reentrant percolation behaviour from 
a non-percolated to a percolated fluid and back to to 
a non-percolated fluid. The addition of 2^ particles to 
the mixture reduces branching, but as the AB bonds are 
energetically favourable, with respect to AA bonds, the 
probability of branching can, under some circumstances, 
increase rather than decrease as the composition of 2^ 
particles increases and re-entrant percolation behaviour 
obtains. 

As the pressure is reduced (panel (02) of Fig. [S]) a 
new two phase region emerges in the form of a closed 
loop bounded above (below) by an upper (lower) critical 
point. It is well inside the percolated region and there- 
fore it is a phase transition between two network fluids 
or percolated states. The new two-phase region is not 
connected to the x — axis, clearly showing that the 
condensation is unrelated to the LV condensation of the 
pure 3a fluid. As the pressure is further reduced, the dis- 
tinct two-phase regions merge and, as previously, there 
is a negative azeotrope (panel (C3) in Fig. [5]). Once the 
azeotrope appears, a new LV branch grows bounded by 
a critical point, wlrich moves towards lower temperature 
and density as the pressure vanishes (see Fig. [S](c)). The 
mixture exhibits 'empty' fluid behaviour. 



Let us focus now on the pT projection of this mixture, 
represented in Fig. [Z](c). It is still a limiting case of type 
I-A mixtures (type I with an azeotrope). Interestingly, 
the critical line leaves the critical point of the pure 3a 
fluid and goes initially to higher temperatures, a very 
unusual situation. 



3- (*AA — 0, f*AB — ^*BB — 1 

The last set of Tx phase diagrams in Fig. [5] (fourth 
row) illustrates the extreme case in which A patches 
do not interact (e^^ = 0). As a result, none of the 
two pure fluids undergoes LV condensation. However, 
the bonding interaction between dissimilar patches {^AB) 
is present and drives condensation at pressures below 
pa 2.92 X IQ-^e/w^. Sightly below (panel (di) of 
Fig. [6]) a closed miscibility gap grows inside the perco- 
lated region. Well below p*^"*") (panels (^2) and (d^) of Fig. 
[5]) there is a negative azeotrope with two LV branches 
bounded by critical points. As the pressure vanishes, 
both critical temperatures tend to zero, but the criti- 
cal density behaves differently in each branch (see Fig. 
[8] (d)). The critical point rich in 2^ particles exhibits 
'empty' fluid behaviour, while the critical point rich in 
3a particles tends asymptotically to r\c — > 0.007. 

The behaviour is analogous to the mixture with = 
0.5 with the absence of the LV branch connected to the 
pure 3a fluid condensation. It is not possible to clas- 
sify this mixture following van Konynenburg and Scotliii 
as neither of the pure fluids has LV condensation. The 
critical line in the pT projection (see Fig. [7](d)) forms 
a closed loop starting and ending at p — s- and T — > 
where we have found two different solutions. 



1 n 



D. 2b - 3a mixtures: (.aa = ess 

Finally, we analyse mixtures with = ^*bb- Rep- 
resentative examples of Tx phase diagrams at p\ = 
pVs/e = 1-05 X 10~^ are plotted in Fig. [HI the pressure- 
temperature projections are depicted in Fig. [TU] and the 
properties of the critical points are shown in Fig. 1111 



■ ^AB > J^) ^AA — ^BB 

In this case AB bonds are energetically favourable with 
respect to both AA or BB bonds, and therefore there is 
no LL demixing at low temperatures. 

The Tx phase diagrams are qualitatively the same as 
those of the reference 2a — mixture. An example is 
plotted in panel (a) of Fig. [9] for t*Ag = 1.1 at p*. It 
should be compared to the phase diagram of the 2 a — ^a 
reference mixture depicted in panel of Fig. [S] At 
constant pressure, the LV two-phase region shrinks as 
the strength of the interaction between dissimilar patches 
increases, and the critical point moves towards higher 
temperature and lower composition. These mixtures are 
limiting cases of type I (Fig. [TOl), and, as expected, the 
critical packing fraction and the critical temperature tend 
to zero as the pressure vanishes (Fig. Ilip . In other words, 
the mixture exhibits 'empty' fluid behaviour. 

Changes in the topology of the phase diagram (as 
the appearance of negative azeotropes) may arise if the 
bonding interaction between dissimilar patches is much 
stronger than the interaction between identical patches: 
(*Agle*aa >> 1, « = A,B. Nevertheless, as far as 
the critical properties at low pressure are concerned, no 
changes are expected to occur. 



2. < e*AB < 1, e*AA = ^"bb 

A small decrease in the AB bonding interaction drives 
phase separation at low temperature, preempting the 
'empty' fluid regime. We have plotted two phase dia- 
grams at the same pressure (p*) in Fig. [5] for: e^s 
0.955 in panel (b) and = in panel (c). For bonding 
strengths < 1 there are two regions of phase sepa- 
ration (LL demixing bounded by an upper critical point 
at low temperatures, and LV condensation bounded by 
a lower critical point at high temperatures) at pressures 
below the critical pressure of the pure 3^ fluid (&). When 
^AB << 1) both regions merge into a single two-phase re- 

(3) 

gion without critical points below pc , see for example 
panel (c). In both cases, at pressures above pc , LL 
demixing is always present bounded by an upper criti- 
cal point. These mixtures are limiting cases of type III 
mixtures (see Fig. [T0|) . The strength of the AB interac- 
tion leads to slight changes of the critical line. Note, for 
example, the slope of the critical line at high pressures: 
Positive when << 1 and negative when e*^g ^ 1. 
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FIG. 9. Phase diagrams of 2s — 3a binary mixtures in the 
scaled-temperature vs composition plane at constant pressure 
Pi = pVs/€ = 1.05 X 10"*. In all cases e^^ = £%b = 1- 
eAB = 1-1 (i)i ^ab = 0.955 (fe), and e^s = (c). The inset 
in (a) is a zoom of the two-phase region. 
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FIG. 10. Pressure-temperature projections of the phase dia- 
grams of binary mixtures with e^^ — egg = 1. 
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FIG. 11. Scaled-critical temperature vs critical packing frac- 
tion of binary mixtures with — e*gg = 1. In all cases with 
^*AB < 1 the critical packing fraction tends to 77c — >■ 1. 



IV. CONCLUSIONS 

We have carried out a systematic analysis of binary 
mixtures of patchy particles with 2 patches of type B 
and 3 patches of type A, focusing on the 'empty' fluid 
behaviour discovered in one of these mixturesii^. We 
have found a rich phase behaviour with closed miscibil- 
ity gaps, liquid-hquid demixing, negative azeotropes and, 
from a topological point of view, three different types of 
mixtures: Type I, I-A and III. The main results are sum- 
marised in Table H] 

We have found that liquid-liquid demixing, at low tem- 
peratures and pressures, is the primary mechanism that 
preempts the 'empty' fluid regime in some of these mix- 
tures. LL demixing is always present when unlike {AB) 
bonds are weaker than at least one of the other two bonds 
between particles of the same species {AA or BB bonds) . 
These mixtures are limiting cases of type III. Within this 



mixture 


empty fluids 


type 


section 


(*BB > 1 


no 


III 


HUB 11 


< e^s < 1 


yes 


I 


HUB 31 


ess = 


no 


I 


HUB 21 


<^AA > 1 


no 


III 


line 11 


< < 1 


yes 


I-A 


line 21 




yes 




line 31 


(■AB > 1 


yes 


I 


nil D 11 


< £as < 1 


no 


III 


nil D 21 


Cab = 


no 


III 


nil D 21 



TABLE I. Summary of the phase behaviour of 2b — Sa mix- 
tures: The bonding energy that characterizes the mixture (the 
other two are set to one), the stability of 'empty' fluids and 
the topological classification of the mixture phase diagram. 



class the detailed phase behaviour depends on the values 
of the bonding energies. For example, the slope of the 
critical line in the pT projection, the appearance of closed 
miscibility gaps or the asymptotic values of the critical 
temperature and critical concentration as the pressure in- 
creases, are controlled by the specific values of the bond- 
ing energies. 

The mixtures exhibit 'empty' fluid behaviour if the AB 
bonds are stronger than at least one of the other two. 
These mixtures are limiting cases of type I or type I-A 
(type 1 with the addition of a negative azeotrope)^. The 
most interesting case occurs in mixtures with bonding 
strengths = and << 1. In these mixtures 
there is a negative azeotrope and the critical line goes 
initially to temperatures above that of the pure 3^1 fluid. 

There is only one case where AB bonds are ener- 
getically favourable and the mixture does not exhibit 
'empty' fluid behaviour: The mixture with ^ab = ^aa 
and ess = 0. In this case the absence of bonding be- 
tween particles with 2 patches precludes the 'empty' liq- 
uid regime. In fact, near the critical point of the mixture 
there is coexistence between two network fluids with fi- 
nite packing fractions as the pressure vanishes. This be- 
haviour is similar to that found previousljsii for 2a — 5a 
binary mixtures {i.e mixtures of particles with 2 and 5 
identical patches) where the entropy of bonding drives 
liquid-liquid demixing preempting the liquid- vapor phase 
transition. 

In a previous studyii we found that the driving force 
for the topological change in the phase diagram of mix- 
tures with a single bonding energy is the difference in 
bonding entropy associated with bonds between like- and 
unlike- particles. It is desirable to address this question 
in more detail by investigating mixtures of particles with 
different diameters, and/or, different types of bonding 
sites in such a manner that both fluids are constrained 
to have similar values of the critical pressure and tem- 
perature. In that context, the present study is a step to- 
wards elucidating the competition between bonding en- 
tropy and bonding energy in the determination of the 
topology of the phase diagrams of binary mixtures. 
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Binary mixtures of particles with different types of 
bonding sites are also crucial to the investigation of a 
wide range of structured fluids, with different macro- 
scopic properties, including stable bigel phases^l. This 
is an immediate goal that will be addressed in future 
work. 

Finally, a word about the absolute stability of the 
fluid phases reported in this paper. Although at low 
pressures and low (but finite) temperatures we expect 
the fluid to be absolutely stable it is clear that solid 
phases will intervene and preempt some of the high pres- 
sure and/or low temperature features of the phase dia- 
grams reported here. The calculation of solid phase dia- 
grams is much more complicated but progress has been 
reported recently, including ground state analysis of bi- 
nary mixtures of low- functionality patchy particles in two 
dimensions^^i^. 



V. ACKNOWLEDGMENTS 

This work has been supported, in part, by the 
Portuguese Foundation for Science and Technology 
(FCT) through Contracts Nos. POCTI/ISFL/2/618 and 
PTDC/FIS/098254/2008, by the R&D Programme of 
Activities (Comunidad de Madrid, Spain) MODELICO- 
CM/S2009ESP-1691, and by the Spanish Ministry of Ed- 
ucation through grant FIS2008-05865-C02-02. D. de las 
Heras is supported by the Spanish Ministry of Education 
through contract No. EX2009-0121. 

ip. N. Pusey, [jTThys . Cond. Matt.'20, 494202 (2008). 

2E. Zaccarelli JJ. Phy s. Cond. Matt. 19, 323101 (2007). 

^M. E. Leunissen, C. G. Christova, A. -P. Hynninen, C. P. Royall, 
A. I. Campbell, A. Imhof, M. Dijkstra, R. van Roij, and A. van 
Blaaderen. [Naturel 437. 235 (2005). 

*A.-P. Hynninen, J. H. J. Thijssen, E. C. M. Vermolen, M. Dijk- 
stra, and A. van Blaaderen, N at Mater,. 6. 202 (2007). 



^A. Stradner, G. FofR, N. Dorsaz, G. Thurston, and P. Schurten- 
berger, Phy s. Rev. Lett.| 99, 19 8103 (2007). 

®A. B. Pawar and I. Kretzschmar, [Macromol. Rapid Commun.| 31, 
150 (2010). 

'^E. Bianchi, J. Largo, P. Tartaglia, E. Zaccarelli, and F. Sciortino, 
JPhys. Rev. Lettr f97, 168301 (2006). 

°E! Bianchi, PT Tartaglia, E. Zaccarelli, and F. Sciortino, 
'J. Chem. Phys. 1 128, 144504 (2 008) . 

''M. S. Wertheim, [j. Stat. Phys^ aS, 35 (1984); 35, 19 (1984); 42, 
459 (1986); 42, 477 (1986). 

lOp. J. Flory rrXm. Chem. Soci 63. 683 (1941); W. H. Stock- 
mayer, J. Chem. Phys. 11, 45 (1943); P. J. Flory, Principles of 
Polymer Chemistry (Cornell University Press, Ithaca, 1981). 

i^D. De las Heras, J. M. Tavares, and M. M. Telo da Gama, Sub- 
mitted to Soft Matter(2010), arXiv: 1012.3605vl 

^ ^J. M. Tavare s, P. I. C. Teixeir a, and M . M. Telo da Gama, 
jPhys. Rev."¥l 80, 021506 (2009); |Mol. Phys 107, 453 (2009). 

^^J. M. Tava res, P. I . C. Teixeira, M. M. T. da Gama, and 
F. Sciortino, [J. Chem. Ph ys. 132, 234502 (2010). 

I' ^J. M. Tavare s, P. I. C. Teixeira, and M. M. Telo da Gama, 
JPhy s. Rev."g |81, 010501R (2010). 

i^JTRusso, J7 M. Tavares, P. I. C. Teixeira, M. M. Telo da 
Gama, and F. Sciortino, Phys. Rev. Lett., accepted for publi- 
cation. (2011). 

^^W. G. Chapman, G. Jackson, and K. E. Gubbins, pol. Phys.| 65, 
1057 (1988). 

^^P. H. Van Konynenburg and R. L. Scott, 

IPhil. Trans. R. Soc. L ond. A 298, 495 ( 1980). 

l^N. F. Carnahan and K. E. Starling, |J. Chem. Phys.| 51, 635 
(1969). 

^^J. Rowlinson and F. L. Swinton, Liquids and Liquid Mixtures 
(Butterworths, England, 1982). 

There is one exception, when = none of the pure flu- 

ids undergoes LV condensation, and therefore it is not possible 
to classify this mixture within the standard classification of van 
Konynenburg and ScotliiX. 

2iA. Goyal, C. K. Hall, and O. D. Velev, |J. Chem. Phys.| 133, 
064511 (2010); Soft Matt er 6, 480 (2010).^ 

^^G. Doppelbauer, E. Bianchi, and G. Kahl, |J. Phys. Cond. Matt.| 
22, 104105 (2010). 

P. K. Doye, A. A. Louis, I.-C. Lin, L. R. Allen, 
E. G. Noya, A. W. W ilber, H. C. Kok, and R. Lyus, 
|Phys. Chem. Chem. Pi^ 9, 2197 (2007). 



